Electronic Entanglement in the Vicinity of a Superconductor 
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A weakly biased normal-metal-superconductor junction is considered as a potential device inject- 
ing entangled pairs of quasi-particles into a normal-metal lead. The two-particle states arise from 
Cooper pairs decaying into the normal lead and are characterized by entangled spin- and orbital 
degrees of freedom. The separation of the entangled quasi-particles is achieved with a fork geom- 
etry and normal leads containing spin- or energy-selective filters. Measuring the current-current 
cross-correlator betw een the two normal leads allows to probe the efficiency of the entanglement 
(|cond-mat/0009193| ). 
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The nonlocal nature of quantum mechanics has been 
demonstrated theoretically Q using entangled pairs of 
particles several decades ago. Recently, potential appli- 
cations of this entanglement have been found in quantum 
cryptography g], in quantum teleportation ||, and in 
quantum computing H. It is thus necessary to search 
for practical ways to produce such pairs given a spe- 
cific interaction between particles. While past experi- 
ments have focused on pairs of photons jq] propagating 
in vacuum, attention is now turning to electronic systems 
||, where this entanglement interaction can be stronger 
while coherence can still be maintained over appreciable 
distances in mesoscopic conductors. A scheme was re- 
cently presented J7J] which discussed the entanglement of 
electrons via the exchange interaction in pairs of quan- 
tum dots. Here, we propose a rather robust electronic 
entanglement scheme based on the Andreev reflection of 
electrons and holes at the boundary between a normal 
metal and a superconductor. 

The basic concept underlying the microscopic descrip- 
tion of superconductivity is the formation of Cooper 
pairs. A normal metal in vicinity to a superconduc- 
tor bears the trace of this phenomenon through the 
presence of Bogoliubov quasi-particles, or through the 
non-vanishing of the Gor'kov Green function | F = 
(ck|C_ki) (cko- denote the usual electron annihilation op- 
erators). While in a superconductor F — A/ A is a con- 
sequence of a nonzero gap parameter A (A is the pair- 
ing potential), the coherence surviving in the adjacent 
normal metal can be understood through the presence 
of evanescent Cooper pairs. These involve two electrons 
with entangled spin- and orbital degrees of freedom, car- 
rying opposite spins in the case of usual s-wave pairing 
and with kinetic energies above and below the supercon- 
ductor chemical potential. This proximity effect has been 
illustrated in several recent experiments H . 

In order to detect this entanglement and implement it 
for applications, it is necessary to achieve a spatial sepa- 



ration between the two constituent electrons. The entan- 
glement apparatus which is proposed here consists of a 
mesoscopic normal-metal-superconductor (NS) junction 
with normal leads arranged in a fork geometry (see Fig. 
1). Using appropriate spin- or energy-selective filters in 
the two normal leads the quasi-particle pairs are prop- 
erly separated and their entanglement can be quantified 
through a comparison of the intra- and inter- lead noise. 
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FIG. 1. Normal-metal-superconductor (NS) junction with 
the normal-metal lead arranged in a fork geometry, (a) With- 
out filters, entangled pairs of quasi-particles (Cooper pairs) 
injected into the lead N3 propagate into leads Ni or N2 ei- 
ther as a whole or one by one. The ferromagnetic filters in 
setup (b) enforce the separation of the entangled spins, while 
the Fabry-Perot type interferometers in the setup (c) separate 
electron- and hole type quasi-particles. 

We start by noting that for a single channel NS wire 
the zero frequency fluctuations of the currents carried by 
electrons with different spins are completely correlated, 

«/<,/-„ ))„=o = [ dt ((I a (t)I_ a (0))) = «J2» W=0) (1) 



hence (((I a — ^-<r) 2 ))w=o = (((•••}} implies the subtrac- 
tion of the average currents) . This perfect correlation in 
the (subgap) motion of the quasi-particles with different 
spins is a consequence of the entanglement of the Cooper 
pairs injected into the normal wire. 

Next, recall that the current noise cross-correlations in 
a SN-NN fork geometry - without filters on the normal 
probes - are positive when the transmission between the 
superconductor and the normal leads is low |10|. The 



unusual sign (for fermions) of these correlations is due to 
paired electrons penetrating the two normal leads sepa- 
rately, c.f., Fig. 1(a). The positive correlations are fur- 
ther enhanced when the competing channel (with paired 
electrons entering the leads jointly) is suppressed through 
the addition of appropriate spin- or energy selective fil- 
ters to the normal leads, see Figs. 1(b) and (c). The 
wave function of the entangled states generated with ideal 
spin/energy filters then is of the type 



« = ak(7;-e, 



■f3\Te±a;±eT<j), (2) 



where the first (second) argument in \</>i;(j>2) refers to 
the quasi-particle state in lead 1 (2) evaluated behind 
the filters and the upper (lower) signs refer to the setup 
projecting the spin (energy); the coefficients a and (3 can 
be tuned by external parameters, e.g., a magnetic field. 
In such a multi-terminal device, the measurement of 
the zero frequency noise cross-correlator {{I a il-<j2 ))w=o 
then serves to detect electron entanglement, in analogy 
with the above single channel NS wire. 

A step like dependence of the gap parameter at the NS 
interface is assumed, subgap transport is specified, while 
the normal leads are single-channel and ballistic. Using 
the scattering formulation of NS transport jll| , the cur- 
rent operator per spin in normal lead n is defined as 
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dede' {[u* ea (x) d x u e > a >(x)] 7e Q 7e'a' 



[v* a (x) d x v e > a >(x)]7 e > a >7l a } exp[i(e - e')t], (3) 



where the operators j ea describe electron and hole Bo- 
goliubov quasiparticles (with positive energies e) on the 
normal side with a = {p, c, n} a multi-index characteriz- 
ing the 'charge' p (= e, h), spin a (= ±1/2), and incidence 

(lead n); / dx 9 = fd x g — gd x f. The associated wave 
functions (u ea (x) and v ea {x)) are expressed in terms of 
the scattering matrix s a . a '] e.g., for an electron with spin 
a incident from lead n and observed in lead m at x m 
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with wave numbers k± = -y/2m(/i s ± e) and the quasi- 
particle velocities v± = hk±/m (p, s is the chemical po- 
tential in the superconductor). The difference between 
the two wave numbers k± will be neglected (/x s 3> A). 



Let us now turn to the fork geometry of Fig. 1: the 
current leads Ni and N2 are connected to N3, which 
itself is terminated with the NS interface. We discuss 
two schemes: a) two ferromagnetic metal contacts (with 
magnetizations in opposite directions) in leads N\ and 
iV 2 block the propagation of the opposite spin, see Fig. 
1(b), b) two energy filters in Ni and N2 (coherent quan- 
tum dots) select the kinetic energies of electrons and 
holes symmetrically above and below the superconduct- 
ing chemical potential (Fig. 1(c)). In both proposals, the 
penetration of a Cooper pair into a given lead is prohib- 
ited, while allowing the split pair to pass the filters. E.g., 
one electron propagates through N\ with spin 'up', while 
simultaneously the other electron (with opposite kinetic 
energy) propagates through N% with spin 'down'. 

The scattering matrix s QjQ ,' has to account for all mul- 
tiple scattering processes: the Andreev reflection at the 
N3S interface can be specified in terms of the trans- 
mission and reflection amplitudes £13, £23, and ru (i = 
1,2) describing the normal-metal part of the device p2| . 
The latter include the scattering by the 'beam splitter' 
N3 < — ► Ni, N2, and account for the presence of spin 
or energy filters in leads Ni and N2. The corresponding 
transmission and reflection amplitudes are found itera- 
tively following the scheme sketched in Fig. 2 and ac- 
counting for all interference processes in the device. 
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FIG. 2. The scattering matrix s a ^ a i determining the noise 
correlator (fUl) incorporates all the internal scattering fea- 
tures of the fork, the beam splitter (— > i (2) , r (2) ), the stubs 
(— » t (1) , r (1) ), and the normal scattering at the NS interface 
f-t, r). 



Beam Splitter: Time reversal invariance is assumed for 
simplicity. It is then possible to express the transmission 
probability, say, between leads Ni and N2 in terms of the 
other two transmissions: 



T$ = T$T™[2 - T^ ± 2(1 - T™)V 2 ]/T™ 7 



(5) 



where T^ ' 
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|£g'|2. For a fully 
4/9. Ref. H focused on 
a setup which is symmetric between land 2. The lower 



symmetric beam splitter T- 



sign in (H) allows to consider also the case where all trans- 
missions are small, thus a more complete parametrization 
is obtained. The reflection amplitudes of the splitter are 



(2) _ ,(2)*,(2).(2) 
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The phases of the reflection and transmission amplitudes, 
in addition to describing the intrinsic properties of the 
links between the three wires, account for the choice of 
the origin in each lead. Practical choices for these origins 
are: a) the position of the NS interface in lead N3 and 
b) that of the filters in Ni and N2 . Typical phases accu- 
mulated during free propagation are exp(±ifc-t&) with b 
a typical length in the beam splitter. 

Stubs: When blocking the propagation through N2 
with an ideal filter, the transmission and reflection am- 
plitudes £^3, r 3 3 , and r^l follow from the transmission 
and reflection amplitudes i^', tf^, £33, r^', and r^f of 
the three bare leads via 



,(1) _ ,(2) _ /,(2) ,(2) _ ,(2) (2)\ //-, 
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(7a) 
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When blocking occurs at N% , the amplitudes are obtained 
by exchanging the lead indices in ( [7a]) and (7b). 

NS boundary: The NS boundary is split into two parts 
describing normal and Andreev scattering (Fig. 2). In or- 
der to include the normal scattering component the two 



scatterers {£ ( 13 ' , ry ,r\{} and {£ NS , r NS , r' NS 
to obtain the next level amplitudes: 



*13 — t 13 t NS /(l 
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} are combined 

(8a) 
(8b) 
(8c) 



where the primed reflection amplitude r^ s is associated 
with the particle incident from the superconductor. Fi- 
nally, we account for perfect Andreev reflection at the NS 
interface; for an electron incident from Ni the scattering 
amplitude reads 



\S e crl,h-cr2\ 



1*13+1 1*23- 



l + kii+| 2 h2-| 2 + 2Re(r 33 +i 



(9) 
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where the indices ± indicate that the energy dependent 
scattering amplitudes have to be evaluated at the pos- 
itive (negative) value of the quasi particle energy e(-s). 
For each energy e only one of the two leads, Ni or N2, 
is open, resulting in a two-terminal device, thus the re- 
lations \ru\ 2 = 1 — |t;3| 2 = I ^33 1 2 hold. The main fea- 
ture contained in Eq. (0) are the Andreev type reso- 
nances building up within the normal-metal leads. These 
resonances are determined through the sign changes in 
Re(r33 + r33_) and their distance ~ hv F /L is determined 
via the Fermi velocity v F and the characteristic size L of 
the region. In addition, zeros appear in the spectral den- 
sity which are a consequence of a vanishing transmission 
for electrons or holes in this three lead geometry. 



The above scheme fully specifies the scattering ma- 
trix s a . a > for the case with ideal filters. For non-ideal 
filters the stub should be replaced with a proper descrip- 
tion of the lead N2 including its non-ideal filter; in ad- 
dition, the normal scatterer described through the am- 
plitudes {*i3, ^3, ri 1} above has to be combined with 
an additional scatterer describing the non-ideal filter in 
the lead Nj. E.g., an energy filter requires inclusion of 
a Fabry-Perot interferometer characterized through the 
scattering amplitudes *i,*2,?*i,?*2 and the separation d 
of the double barrier system and producing a transmis- 
sion t FP — £it 2 exp(ifc<i)/[l — rir' 2 exp(2ikd)]. The reso- 
nance spacing should be larger than the applied bias for a 
proper device operation as a filter. The initial resonance 
lines produced by the quantum dot will then be deco- 
rated by Andreev-type resonances and zeros originating 
from the NS-fork structure. 

The above entangler is essentially a two terminal de- 
vice where electrons with, say, a given spin from lead 1 
are converted into holes with an opposite spin in lead 
2. The current correlations between 1 and 2 are positive 
and can be obtained using the definition of the noise in 
combination with (|3|), at T = 0, 

«W-ri)) = j5]/ l de\s a , a ,\ 2 (l-\s a , a ,\ 2 ), (10) 



where a voltage eV\ is applied between the lead Ni and 
the superconductor while keeping the lead N2 unbiased. 
For the case of ferromagnetic filters, the chemical poten- 
tial which enters the (sharp) electron and hole distribu- 
tion functions depends also on the spin index. The multi- 
indices a and a' to be summed over in (|1C| ) depend on the 
type of filters in the normal leads Ni and N2: For ferro- 
magnetic filters (SN-FF) with the spin in F X ( 2 ) pointing 
up (down) a = {e(h) | 1} and a' = {h(e) J. 2} (the prop- 
agation of other states is blocked by the filters). On the 
other hand, for the setup selecting a definite quasi par- 
ticle energy via Fabry-Perot type filters we have to sum 
over spins with a — {e 1 (J.) 1} and a' = {h I (|) 2} (we 
assume filters selecting quasi particles and quasi holes 
in leads Ni and N2, respectively). Applying the same 
voltage to the lead N2 as well does not change the an- 
swer in the normal fork (SN-NN) but renders the result 
for the ferromagnetic filters (SN-FF) twice larger. The 
current fluctuations (UQ) are straightforwardly converted 
into "counting" correlations as known from quantum op- 
tics: with en an {t) = L dt 1 I an (t') we find: 



((n CT i(*)n_ CT2 (t)»|t 



(i/e 2 )«/ CT i/- CT 2)}c=o . (11) 



and hence ({{n al - n_ . 2 ) 2 ))/((" 2 ri)) ~ 0. The result 
( p0| ) together with the fact that the two currents I^i and 
I-u2 are necessarily correlated constitutes the main jus- 
tification for the proposed entanglement device: Eq. (nfl) 



corresponds precisely to the current noise in lead 1. Ide- 
ally, these correlation measurements should be performed 
using fast electronics in order to generate time resolved 
voltage pulses for electron injection/detection. 

This electronic entanglement apparatus can now be 
completed with a detection apparatus in order to test 
non-local correlations (Bell inequalities). For an entan- 
gler based on ferro/energy filters, the detection apparatus 
involves filters of opposite type (energy/ferro) . Concen- 
trating on energy filters, a positive energy particle emerg- 
ing in lead Ni can have either spin orientation, which can 
be "measured" by connecting this lead to, say, a magnetic 
contact with known spin orientation. In the opposite arm 
one should have a similar contact but with a magnetiza- 
tion axis rotated by n/2 in order to achieve the analog of 
the spin correlation experiments of Ref. [ 5L A detailed 
discussion of these will be provided later [ 1J| . 

The proximity induced entanglement of quasi particles 
in NS-fork type devices was implicit in Refs. p5[ and Ref. 
[|§. Consider the above SN-NN setup with the lead Ni 
biased with respect to the superconductor, while keeping 
the lead N2 at the superconductor chemical potential. 
A finite current IziYi) = (2e//i) j Q 1 de |s e fi,/ij.2| will 
flow through lead N2 in response to the bias eV\ across 
lead Ni . While both experiments in Ref. JL5| use a mag- 
netic field to separate electron- and hole type quasi par- 
ticles, the more recent suggestion in Ref. Q proposes 
two ferromagnetic needles, a setup similar to our SN- 
FF device. This Andreev drag effect is quite robust and 
decreases only gradually with decreasing quality of the 
filters. The condition for this drag effect to persist reads: 
J e 1 de[\s enthi2 \ - |seTi,eT2| ] > 0, implying that the 
normal current injected from lead Ni to lead N2 remains 
smaller than the 'drag current' due to Andreev reflected 
holes. Note that the current I2 will vanish when replac- 
ing the s-wave superconductor with a p-wave material or 
a normal conductor. 

Summarizing, we have proposed an electronic entan- 
glement device based on the proximity effect and have 
shown how to probe the resulting non-local electronic 
correlations in an emphatic way through a measurement 
of the current cross-correlator. Using a special fork ge- 
ometry with, say, Fabry-Perot filters one arrives at a 
natural source of spin-entangled electron pairs, a device 
with potential applications in quantum computing archi- 
tectures based on spintronics |l7[ |. This device presents 
the advantage - as compared to its ferromagnetic cousin 
- that it can be realized with nowadays splitters JlS| ] 
and quantum dot technology, e.g., using semiconductor- 
superconductor heterostructures |1SJ], and does not re- 
quire interdot coupling. Moreover, this SN-NN device 
appears to be more promising regarding potential appli- 
cations for quantum information processing: the inser- 
tion of Fabry-Perot filters destroys only the orbital en- 



tanglement of the electrons, while the (most valuable) 
spin entanglement persists, contrary to the situation in 
the SN-FF device where the filters project the spin, but 
where the entanglement of energy degrees of freedom per- 
sists nevertheless. 
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